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Abstract 

In the theory of weakly non-hnear elasticity, Hamilton et al. [J. 
Acoust. Soc. Am. 116 (2004) 41] identified W = iil2 + {A/'i)h + Dll 
as the fourth-order expansion of the strain-energy density for incom- 
pressible isotropic solids. Subsequently, much effort focused on the- 
oretical and experimental developments linked to this expression in 
order to inform the modeling of gels and soft biological tissues. How- 
ever, while many soft tissues can be treated as incompressible, they 
are not in general isotropic, and their anisotropy is associated with the 
presence of oriented collagen fiber bundles. Here the expansion of W 
is carried up to fourth-order in the case where there exists one family 
of parallel fibers in the tissue. The results are then applied to acous- 
toelasticity, with a view to determining the second- and third-order 
nonlinear constants by employing small-amplitude transverse waves 
propagating in a deformed soft tissue. 
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1 INTRODUCTION 



Due to their high fluid content, many soft biological tissues and gels exhibit 
nearly incompressible behavior: they are constrained to undergo essentially 
volume-preserving deformations and motions. This has long been known and, 
in fact, can traced back to the experiments of Jan Swammerdam[Tl[2] in 1758. 
Mathematically, the internal constraint of incompressibility is expressed by 
imposing the condition det = 1 at all times, where F is the deformation 
gradient. In exact nonlinear isotropic elasticity theory, the strain-energy 
density W is written in fullest generality as a function of /, //, ///, the first 
three principal invariants of the right Cauchy-Green strain tensor C = F'^F, 
which are defined by 

/ = trC, II = llitT -tT{C% /// = detC = (detF)2. (1) 

It is then easy to enforce incompressibility, which requires that III = 1 at 
all times and hence that W = W{I, II) only. 

In weakly nonlinear isotropic elasticity theory, W is expanded in terms of 
invariants of the Green-Lagrange strain tensor E = |(G — /), either using 
the set 

ii = tr E, i2 = i[(tr Ef - tiiE^)], i^ = det E, (2) 

(this is the choice of Murnaghan[3]) or the set 

Ji = tr E, h = tr E^, h = tr E^, (3) 

(this is the choice of Landau and Lifshitzjl]). In each case the respective 
terms are of orders 1, 2 and 3 in E. 

Then, in third-order elasticity theory, W is expanded as a linear com- 
bination of ii, 12, is, iii2, and if (or, equivalently, /i, I2, I3, hh, and If) 
and higher-order terms are neglected — and thus, 5 elastic constants are re- 
quired to describe third-order solids. [Note that some authors refer to this as 
second-order elasticity since the stress is second order in the strain; a similar 
comment applies to fourth-order versus third-order elasticity.] In fourth-order 
elasticity theory, terms linear in if, i\, iii^, and 2^22 (or, equivalently. If, J|, 
/1/3, and I1I2) are also retained — and thus, 9 elastic constants are required 
to describe fourth-order solids. Now, however, it is slightly more complicated 
to enforce incompressibility, because the constraint involves a combination of 



2 



invariants of E, rather than a single invariant of C. Nevertheless, this can be 
achieved because of obvious connections between the sets ([T]), ([2]), and (jSJ. 
This was done long ago by Rivlin and Saunders[5] in 1951, and by Ogden[6] 
in 1974, for third- and fourth-order incompressible solids, respectively. In 
particular, these authors showed that the number of elastic constants drops 
from 5 to 2 for third-order solids, and from 9 to 3 for fourth-order solids 
compared with the compressible case. 

Recently, Hamilton et al.[7], unaware of these previous contributions, 
rediscovered the latter result, and wrote W in the form 

W = i^h + ^/s + Dli, (4) 

where /i. A, and D are the respective second, third-, and fourth-order elastic- 
ity constants of weakly nonlinear incompressible isotropic elasticity. Subse- 
quently, several papers p[9l[ia[T3l[Ta[ini[IIl[15l[l6l[I71[l8l[2ll[ia 
have been devoted to the study of the theoretical and experimental impli- 
cations of this expansion for linearized and nonlinear waves and vibrations 
in incompressible solids. Experiments have been conducted on isotropic gels 
and phantoms to measure the nonlinear elastic coefficients, and good progress 
has been achieved in the understanding and implementation of elastographic 
techniques [121 [131 El [HI CI] • 

Now, the ultimate goal of elastography is the imaging of biological soft 
tissues, and it is therefore important to be able to evaluate /i. A, and D 
experimentally. At the same time, it must be kept in mind that soft tissues 
are anisotropic solids, in contrast to gels and phantoms. This anisotropy 
is due to the presence of oriented collagen fiber bundles, which are three 
orders of magnitude stiffer than the surrounding tissue in which they are 
embedded (this surrounding tissue consists of a network of elastin fibers, 
smooth muscle cells and proteoglycans, inter alia, which collectively form an 
essentially isotropic 'matrix'.) The natural conclusion of this consideration is 
that there is a need for expansions of W up to the third or, depending on the 
context, fourth order for incompressible anisotropic solids. For the particular 
application considered in Section III of the present paper, however, it suffices 
to pursue the analysis only up to the third order. 
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2 FOURTH-ORDER TRANSVERSELY 
ISOTROPIC SOLIDS 



Here we consider soft tissues with one preferred direction, associated with a 
family of parallel fibers of collagen. We denote by A the unit vector in that 
direction when the solid is unloaded and at rest. The theory of Spencer [21] 
tells us that in all generality, is a function of at most 5 invariants in exact 
nonlinear compressible elasticity theory. One valid choice for these quantities 
is 



W = W{h,l2,h,h,h), where h = A-EA, h = A-E^A. (5) 



Now we expand W up to terms of order 4 in the Green-Lagrange strain. This 
involves a linear combination of the following quantities: 



Note that the first-order terms are omitted from the list because they give rise 
to constant stresses. Hence there are 29 elastic constants for fourth-order, 
transversely isotropic, compressible solids. 

In incompressible solids, however, the isotropic invariants are not inde- 
pendent. They are connected by the incompressibility condition, which can 
be written in terms of the invariants ([2]) exactly as[n] ii = —2^2 — 4^3, or 
equivalently, in terms of the invariants ([3]), as 



so that the number of independent quantities in ([6]) is greatly reduced. For 
instance. If is of fourth order, /1/2 = /| and If = at fourth-order. This 
allows the list to be shortened to 



order 2: if, I2, I^, h, hh, 

order 3: Jf, hh, h, IJ^ hh, Iih, hh, /|, hh, 
order 4: if, P,h, hh, ll Ul ifh, ll ll ifh. 



(6) 



Il=l2- Ih - ll + 2/1/2 - §/?, 



(7) 



order 2: I2, ij, I5, 

order 3: I3, hh, /|, hh, 

order 4: ll hll hh, ll ll hh, 



(8) 



which are associated with 13 elastic constants in toto. 
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For third-order, transversely isotropic, incompressible, nonlinear elastic 
solids, for example, the strain-energy density is thus expressible in the form 

W = iJ,l2 + —I3 + ctill + "2/5 + ashh + OiJl + ctbhh, (9) 

where /i, ai, 0^2 are second-order elastic constants and A, a^, a^, a^, 
are third-order elastic constants. A lengthier expansion is required at fourth 
order, which requires a linear combination of the six additional terms listed 
in dH]). We do not include these here. Indeed, it suffices for our present 
purposes to restrict attention to second- and third-order terms in W since 
at this order the speed of infinitesimal plane waves depends linearly on the 
strain through the third-order constants, even for the specialization to an 
isotropic material. 



3 ACOUSTOELASTICITY 

We now consider the propagation of small- amplitude plane body waves in a 
deformed soft tissue with one family of parallel fibers. We assume that the 
solid has been pre-stressed by the application of the Cauchy stress cr, say, 
giving rise to a pre-deformation with corresponding deformation gradient F. 

Let u = u{x,t) denote the mechanical displacement associated with the 
wave motion, p the (constant) mass density, and p the incremental Lagrange 
multiplier due to the incompressibility constraint. Then the incremental 
equations of motion and of incompressibility read[ 



•AopigjUj^pg - p^i = pd'^Ui/dt'^, Ui^i = 0, (10) 

respectively, where the comma signifies partial differentiation with respect 
to the current coordinates x = Xi and .Aq is the fourth-order tensor of 
instantaneous moduli, with components 

dW d^w 

•^Opiqj = FpaFqjB^ij—^ h FpaFqj^FjyFi-^-— — ; . (ll) 

This expression is exact, while to compute the derivatives of W in respect of 
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the following quantities are needed: 
dh 



'2 



OEa^dEpy 
dh 



2Eai3, 



dEap 

'3 



dEa^dEpy 2 



{^apEy^ + ^auEp^ + 6^uEa/B + d^jiEau) 



d'h 3 



dEap 

0, 



dEa^dEpy 
dh 



dEap 

'5 



dEa^dEpy 2 



Aqi^A^E -\- AjjA^Eq.^, 

{AaAydp^ + AaApd-^y + A/jA^^ci, + A^Aydap)- (12) 



a^J, 1 



Let us now consider the propagation of homogeneous plane body waves, 
in the form 

u = ae^'^^"-^-^*), p = iA;Pe^'=("-"'-'^*), (13) 

where a is the unit vector in the direction of linear polarization, k is the 
wave number, v is the phase speed, and P is a scalar. Then the incremental 
incompressibility condition f llOp 9 gives: a n = 0, and the wave is thus 
transverse. The incremental equations of motion f lTQ|) i can be written in the 
form 

Q{n)a — Pn = pv'^a, where [Q{n)]ij = AopigjUpUg. (14) 

Taking the dot product of this equation with n gives an expression for P. 
Substituting this back into (fT4|) and using the orthogonality of the propa- 
gation and polarization vectors, we end up with the symmetric eigenvalue 
problem, 

[I — n]Q{n)[I — n n\a = pv'^a (15) 

for the wave speed and polarization for any given direction of propagation. 
The wave speed v is given simply by pv^ = AopigjUpUgaiaj. 
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As a simple example of the application of the above equations, consider 
the case where a sample of soft tissue is under uniaxial tension or compression 
with the direction of tension parallel to the fibers. We denote by e the 
elongation of the sample in the direction [xi, say) of the uniaxial stress and 
discard all terms equal to or higher than in the expansions. In particular, 
we obtain the approximations 

F = diag[l + e, 1 - e/2, 1 - e/2] , E = diag[e, -e/2, -e/2] . (16) 

Also, we find that there are only 15 non-zero components of ^Aq in the 
{xi,X2, X3) coordinate system because, in this special case, the principal axes 
of stress and strain are aligned. Now let 6 be the angle between the direc- 
tion of propagation and the Xi axis. Then the secular equation, which gives 
the wave speed in terms of the elastic moduli, has the same form as in the 
isotropic case [26], with the difference that the moduli depend not only on 
the isotropic elastic constants /i and A, but also on the anisotropic elastic 
constants ai, a^-, as, 04, and 05. It reads [2^ 

p^;2 = (a + ^ _ 2/3) cos^ e + 2(/3 - 7) cos^ ^ + 7, (17) 

where 

01 = "^01212; 2/3 = .4,01111 + ^02222 " 2^01122 " 2^01221; 7 = -4.02121- 

(18) 

Lengthy, but straightforward, calculations reveal that 

2/3 - a - 7 = 2ai + 2(4ai + 3a2 + 3a3 + 804 + 2^5)6, (19) 
2(/3 — 7) = 2ai + (3/i + lOai + 802 + 603 + 604 + 4a5)e, 
47 = 4/i + 2^2 + {A + 2a2 + 4^3 + 2^5)6. 

Clearly, Eq. (fTTl) provides a direct means of evaluating the elastic constants, 
by measuring the transverse wave speed for a variety of e and 6 combinations. 



4 CONCLUSION 

We have shown that 7 elastic constants are required to describe third-order 
incompressible solids with one family of parallel fibers, and we have indicated 
that these constants may be determined by using the acoustoelastic effect. 
We have also noted that 13 constants are needed for fourth-order solids. 
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In closing this Letter, we remark that although many biological soft tis- 
sues exhibit transverse isotropy (one family of parallel fibers), many more 
must be modeled as orthotropic materials (two families of parallel fibers). 
Then, the list of invariants is increased from 4 (/i, /2, /4, I5) to 7, with the 
following 3 additions: 

h = B- EB, h = B E^B, h = {A- EB)A ■ B, (20) 

where B is the unit vector in the direction of the second family of parallel 
fibers. In that case, an expansion of W up to third- or fourth-order yields 
such a large number of elastic constants to be determined, that it defeats 
its own usefulness. We would argue that it is more advantageous to turn to 
the exact nonlinear theory, where constitutive anisotropic models have been 
successfully evaluated against experimental test data to capture well the 
mechanics of orthotropic soft tissues. For instance, the Holzapfel-Gasser- 
Ogden model for arteries [27] is now implemented in many Finite Element 
Analysis software packages. It requires the experimental determination of 
only 3 elastic constants per layer of arterial wall. 
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